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ABSTRACT 

We present the mass-X-ray observable scaling relationships for clusters of 
galaxies using the XMM-Newton cluster catalog of Snowden et ah. Our results 
are roughly consistent with previous observational and theoretical work, with one 
major exception. We find 2-3 times the scatter around the best fit mass scaling 
relationships as expected from cluster simulations or seen in other observational 
studies. We suggest that this is a consequence of using hydrostatic mass, as 
opposed to virial mass, and is due to the explicit dependence of the hydrostatic 
mass on the gradients of the temperature and gas density profiles. We find a 
larger range of slope in the cluster temperature profiles at r 50 o than previous 
observational studies. Additionally, we find only a weak dependence of the gas 
mass fraction on cluster mass, consistent with a constant. Our average gas mass 
fraction results also argue for a closer study of the systematic errors due to 
instrumental calibration and modeling method variations between analyses. We 
suggest that a more careful study of the differences between various observational 
results and with cluster simulations is needed to understand sources of bias and 
scatter in cosmological studies of galaxy clusters. 

Subject headings: cosmology: observations — galaxies: clusters: general — X- 
rays: galaxies: clusters 
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1. Introduction 


Studies of clusters of galaxies provide for a variety of cosmological tests (see Voit 2005, 
for a recent review). The precision of these tests is limited by the accuracy and precision of 
the scaling relations used to transform X-ray observables, such as temperature or luminosity, 
into cluster mass measurements (e.g., Rapetti et al. 2008; Vikhlinin et al. 2009). Improved 
modeling of the complex physics at work in clusters has lead to better study of the sources of 
bias and scatter in the mass scaling relationships. These models suggest that the hydrostatic 
equilibrium assumption used to calculate total cluster masses from X-ray data can lead 
to underestimates on the order of 10-20% (e.g., Nagai et al. 2007b). The discrepancy is 
explained by the presence of nonthermal pressure support in clusters. Theoretical studies 
have found that the expected scatter in mass scaling relationships is dependent on the 
dynamical state of the clusters under study (e.g., Kravtsov et al. 2006). Relaxed systems 
should show lower scatter around the best-fit relation than disturbed clusters. 

This information can be used to tailor observational programs, particularly by limiting 
the bias and scatter in the mass scaling relationships. One suggestion is to use only relaxed 
clusters in cosmological studies. However for high redshift samples, relaxed clusters make up 
only a minority of the cluster population (see e.g., Vikhlinin et al. 2009). Alternatively, the 
right choice of mass proxy could provide low scatter data. Kravtsov et al. (2006) proposed a 
new X-ray proxy for cluster mass, the parameter Y x which is the product of X-ray temper- 
ature and gas mass. The Y x parameter is related to the total thermal energy of the cluster 
gas and is an X-ray analogue to the Sunyaev-Zel’dovich (SZ) flux. Prom cluster simulations, 
Kravtsov et al. (2006) found that the scatter in the cluster mass-lV scaling law is not only 
lower than for other commonly used mass proxies, but shows little dependence on cluster 
dynamical state. 

Observational studies of the Y x parameter as a mass proxy have been limited in the 
number of clusters under study and the statistical quality of the data (Arnaud et al. 2007; 
Vikhlinin et al. 2009). The statistical error has been on the order of the measured scat- 
ter, making it difficult to determine the intrinsic scatter in the relationship. Additionally, 
previous observational work has focused solely on a limited number of relaxed systems. A 
larger and more sensitive observational study of clusters would better test the Y x mass proxy 
and the results of cluster simulations. In this Letter, we use the data from a recent high 
signal-to-noise XMM-Newton survey of nearby clusters to test the usefulness of various mass 
proxies. We assume H 0 = 70 km s _1 Mpc -1 , Qm = 0.3, and 0 A = 0.7. 



-3- 


2. Data Analysis 

Our cluster sample was taken from Snowden et al. (2008), which presented the projected 
temperature, abundance, and surface brightness profiles for 70 clusters found by fitting 
XMM-Newton data (see that work for details of the spectral analysis). While not an unbiased 
sample, the large number of clusters does provide a wide range of cluster properties to study. 
It is not limited to relaxed and/or hot (> 5 keV) clusters. Snowden et al. (2008) discuss 
some of the selection criteria used to produce the sample. We note that highly asymmetric 
clusters and those with strong substructure were excluded. 

We used the projected profiles to determine the gas density and three-dimensional tem- 
perature profiles following the procedure of Vikhlinin et al. (2006). In addition, we used 
the same gas density and temperature models with some simplification. The lower spatial 
resolution of our data did not constrain the second /3-model component required to fit the 
Chandra data in the cluster center, therefore we did not include it. In addition, we compared 
fits with and without steepening at large radii and excluded the steepening when it did not 
produce a significantly better fit of the data. For the temperature fits, when cooling in the 
core was not obvious we fixed the cool component parameters to produce a flat profile in the 
center, i.e. T min /T () = 1 and a coo i = 0. 

We calculated the total cluster mass distribution for each cluster assuming hydrostatic 
equilibrium and the radii (r 2 5 oo and r 500 ) where the cluster mass equaled 2500 and 500 times 
the critical density at the cluster redshift. We do not include clusters where the calculated 
r 25 oo and r 500 values extend beyond the radial coverage of our data. We find that 60/70 
clusters have data extending to at least r 2 500 and 28/70 have data out to r^oo ■ We determined 
the gas mass ( M g ) and total cluster mass (M) enclosed by r 250 o and r 50 o for each cluster. 
We also calculated the average spectral temperature (Tx) within the 0.15-1 r 5 00 radial range 
using the formulation of Vikhlinin (2006). 

Uncertainty intervals were obtained from Monte Carlo simulations. We simulated sur- 
face brightness and projected temperature distributions by scattering the observed data 
according to the measurement uncertainties found in Snowden et al. (2008). These data 
were then fit with the gas density and temperature models and a full analysis performed to 
determine M, M g , and Tx- The uncertainties on these were obtained from their distribu- 
tion in the simulated data. For values evaluated at r 50 0 (r 250 o), the uncertainty includes the 
uncertainty on r 500 (r 25 oo)- 



-4 - 

3. Comparison of Scaling Relations 

We determined the best-fit scaling relationships at r^o using E(z) n M = C(X/X 0 ) Q for 
X = Tx, M g , and Yx- We fixed n — 1, 0 and 2/5, as consistent with our cosmology, and 
X 0 = 5 keV, 4 x 10 13 M 0 and 3 x 10 14 M 0 keV for the Tx, M g and Yx fits, respectively. We 
also present the best- fit relationship between the gas mass fraction, f g = M g /M, and M, as 
characterized by the equation f g = / 3i0 -P a\og w [M / \Q lb M & ] (see Vikhlinin et al. 2009). 

We use the BCES fitting routines 1 which provide a linear regression algorithm that 
allows for intrinsic scatter and nonuniform measurement errors in both variables (Akritas 
& Bershady 1996). We find that the Y|X and orthogonal slope estimators provide the 
most significant (and consistent) results compared to other methods, including bisector. We 
therefore only present the results of the Y|X and orthogonal methods for each of our fits (see 
Table 1). We also include the best-fit results when the relationship slope, a, is fixed at the 
expected value from self-similarity. 

For the M — X relationships we estimate the intrinsic scatter using a generalized form 
of the estimated scatter, 8M/M, used by Vikhlinin et al. (2006): 

/ 5M \ 2 1 y. [Mj - CiXi/XoTl - A M[ 

V M j JV-2^ Mf ’ U 

where A Mi are the measurement errors. Similarly for the f g — M relationship, we calculate 
the scatter 8f g /f g by: 

(Sf t V 1 y— [/g,i — (/g,0 + Q logio[:y/10 15 Mol)] 2 — A/g,i , 

\fj N-2^ fl ' U 

where A /^, are the measurement errors. To compare our scatter with logarithmic scatter 
estimates (e.g., Jeltema et al. 2008; Pratt et al. 2008), multiply the logarithmic estimates by 
In 10 = 2.30. Given the high statistical quality of our data, the scatter in the relationships 
is dominated by intrinsic scatter. 


3.1. Mass— Temperature 

Table 1 and Figure 1 show the best-fit parameters for the M — Tx scaling relation. Our 
best fit (log 10 C* = 14.64, a = 1.67) is consistent with the scaling relations found in other 
datasets (Arnaud et al. 2007: Vikhlinin et al. 2009, log 10 C = 14.580 and 14.635, and a = 


Routines available at http://www.astro.wisc.edu/~mab/archive/stats/stats.html 
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1.71 and 1.53, respectively). We note that Arnaud et al. (2007) use a different definition of 
Tx which integrates the observed temperature profile over 0.15-0.75 r^oo- Since the average 
cluster temperature profile falls at large radii (see e.g., Leccardi & Molendi 2008, Juett et al. 
2009, in prep), this definition will produce higher values of Tx and subsequently lower values 
of the M — Tx normalization. We find that (Xx(0.15 — 0.75r 5 oo)/Tx(0.15 — lrsoo)) = 1-06, 
which translates to a reduction of log 10 C of 0.04-0.05, enough to explain the discrepancy 
between our results and Arnaud et al. (2007). 

When compared to theoretical calculations (Kravtsov et al. 2006; Nagai et al. 2007a, b; 
Jeltema et al. 2008), our results are consistent when the hydrostatic mass is considered, 
but our normalization is lower by ~ 20% when compared to scaling results that use the 
true cluster mass (14.70-14.75). This is a well known result and is likely due to non-thermal 
pressure support that is not accounted for in the hydrostatic mass estimate (e.g., Nagai et al. 
2007a). The biggest difference between our results and previous studies is the difference in 
scatter around the best-fit scaling relation. We find a scatter of 43% which, due to our 
systematics limited dataset, can be attributed to the intrinsic scatter of our sample. Other 
studies, of both observed and theoretical cluster samples, have found significantly lower 
values for the intrinsic scatter, « 10 — 25% (Vikhlinin et al. 2006; Arnaud et al. 2007; Nagai 
et al. 2007a, b; Jeltema et al. 2008). 


3.2. Mass— Gas Mass 

The gas mass-total mass scaling relationship has been previously characterized in two 
ways: a powerlaw scaling between M and M g , and more recently by a linear scaling of f g 
and the logarithm of M. 

Our M — M g results (log 10 C = 14.503 ± 0.018), are in marginal agreement with the 
normalization found by Arnaud et al. (2007, 14.542±0.015). Our best-fit slope is also steeper 
(0.93±0.05 versus 0.80±0.04). Again we find that theoretical models that use the true cluster 
mass have higher predicted normalizations but our results are consistent when hydrostatic 
mass estimates are used (Kravtsov et al. 2006; Nagai et al. 2007a). The slope estimate 
is consistent with theoretical results for both true and hydrostatic mass estimates. Our 
15% scatter is close to the ~10% scatter found in both observational and theoretical work. 
Interestingly, our results match the combination of normalization, slope, and scatter found in 
the simulations of Nagai et al. (2007a) when hydrostatic mass is used, however we differ from 
their Chandra observational results. Our slope is larger (0.93 vs 0.81) while our normalization 
is lower (14.503 vs 14.59). This discrepancy may be due to differences in the calibration of 
the instruments that has been previously noted (e.g., Snowden et al. 2008). 
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We find that the normalization of the f g -M relationship, 0.134±0.005, is consistent with 
the value of Vikhlinin et al. (2009, 0.130±0.007). Our data are consistent with a constant 
slope over the range of mass considered, while the Vikhlinin et al. (2009) result prefers a 
reduction in f g for lower mass clusters. We note that the Vikhlinin et al. (2009) result is 
consistent with gas mass fraction results from groups of galaxies (Sun et al. 2009). Arnaud 
et al. (2007) suggested that the gas mass fraction may be constant above 2 - 3 x 10 14 M 0 , 
and then dropping at lower masses. This would explain both our result and the lower gas 
mass fractions seen in groups of galaxies. 

The difference in mean gas mass fraction between the Vikhlinin et al. (2006) work, 
a subset of the Vikhlinin et al. (2009) sample, and ours is 10-20%. At r 50 o, we find a 
mean f g = 0.1323 ± 0.0019, compared to 0.110T0.002 for the Vikhlinin et al. (2006) sample 
(adjusted to account for differences in cosmology). If we restrict our analysis to clusters with 
kT > 5 keV, the difference is <10%, 0.138T0.003 for our work compared to 0.123T0.003 
(Vikhlinin et al. 2006). 

Our results are also close to those found by Allen et al. (2008). At r 250 o, Allen et al. 
(2008) found a mean cluster gas mass fraction of 0.113T0.003 for clusters with kT > 5 keV 
and low redshifts (z < 0.15). For their full sample (kT > 5 keV but all redshifts) they find 
a mean f g = 0.1104 ± 0.0016. Vikhlinin et al. (2006) noted that their mean f g at r 2 500 was 
significantly less (0.091±0.002, a ~25% difference) than an earlier (but consistent) Allen et 
al. sample. We find fg^oo = 0.1057 ± 0.0005 for clusters with kT > 5 keV, a 4% difference 
with the Allen et al. (2008) results and a wl5% difference with the Vikhlinin et al. (2006) 
work. 


3.3. Mass— YV 

Kravtsov et al. (2006) suggested that a lower scatter (< 10%) proxy for cluster mass 
is the parameter Y x = T x M g . Table 1 gives our best-fits for the M - Y x scaling relation. 
Our best-fit normalization (14.657±0.018) is consistent with the results of Arnaud et al. 
(2007, 14.653±0.015 when renormalized) and Vikhlinin et al. (2009, 14.684±0.015). When 
comparing with theoretical results, we find consistency with hydrostatic mass results (14.645) 
but not true cluster mass (14.712; e.g. Nagai et al. 2007a). The slope of the best-fit M - Y x 
compares well with both observational and theoretical studies. 

The largest difference comes in the measured scatter in the M - Y x relation. We find 
a scatter of « 22%, much larger than the < 10% level expected from simulations. While 
other observational studies have not found as large a scatter (e.g., Arnaud et al. 2007), we 
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note that ours is the first study to be systematics limited and includes twice the number of 
objects. Thus we are able to measure the scatter without a significant contribution from the 
statistical error. We point out that in theoretical work, scatter increases when hydrostatic 
masses are used compared to true masses, though not as large as our result (Nagai et al. 
2007a; Jeltema et al. 2008, 8-20%). 


4. Correlation of Scatter and Deviation with Cluster Properties 

Scatter in the mass-X-ray observable relationships is an important contributor to the 
total error budget in cosmological studies of clusters (see e.g. Vikhlinin et al. 2009). Given 
our large sample of clusters, we can study what factors are most important in producing 
scatter in these relationships which can then be used to refine cosmological studies to reduce 
the scatter. 

One possible cause of the observed scatter is the dynamical state of the cluster. Relaxed 
clusters are expected to better follow the assumptions of hydrostatic equilibrium. Disturbed 
clusters may have additional pressure and energy inputs due to the merging events and the 
additional complication of asymmetric geometries (see e.g. Nagai et al. 2007b). 

We looked for correlation between the deviation and scatter of the calculated hydro- 
static mass, Mi , from the expected mass given the best-fit scaling relationship, M (X,;), with 
measures of the ellipticity and asymmetry in the cluster images (see Figure 2). We identify 
the deviation as SM(X)/M = [Mi — M(Xi)]/Mi and the scatter as ( SM(X)/M ) 2 . The ellip- 
ticity and asymmetry were calculated following the work of Hashimoto et al. (2007). Relaxed 
clusters are expected to have low ellipticity and asymmetry values, while disturbed systems 
will have higher ellipticity and/or asymmetry values. Hashimoto et al. (2007) showed that 
ellipticity is correlated with the P2/P0 power ratio and that asymmetry is related to the 
P3/P0 power ratio (see e.g., Buote & Tsai 1995; Jeltema et al. 2008, for a discussion of 
power ratios). We find no correlation between the amount of deviation or scatter and either 
ellipticity or asymmetry for any of our scaling laws. 

We looked for other cluster properties that might influence the scatter. Vikhlinin 
et al. (2006) noted that implicit in the calculation of hydrostatic mass, there is a depen- 
dence of the normalization of the M — T x relation on the sum of the temperature gradient, 
Pt = (— 1/3) dlogT/dlogr, and gas density gradient, p e ff = (— 1/3) dlogp/dlogr (see their 
Appendix A). If the spread in values is large enough, this dependence should cause a pre- 
dictable deviation around the best-fit relationship. We find a strong correlation between 
SM(X)/M and /3 e // + Pt for all scaling relationships (Figure 3). The Spearman rank cor- 



- 8 - 


relation coefficients were 0.62, 0.53, and 0.61 for X = T x , M g , and Yx, respectively. No 
correlation was found for scatter and l3 e f / + fit- 

The cluster sample used by Vikhlinin et al. (2006) has a narrow distribution of fi e jf + fi t 
values. Our sample however, shows a large variation in the temperature gradient. For fi t , 
we find a mean of 0.53 and a standard deviation of 0.50. The range of temperature gradient 
values is a reflection of the variation of temperature profiles at large radii in our cluster 
sample (Snowden et al. 2008). The gas density gradient has a narrower distribution with a 
mean fi e ff = 0.63 and a standard deviation of 0.10. Our gas density gradient results are in 
good agreement with those found in the REXCESS study (fi ne (0.3-0. 8 r 500 ) = 0.60T0. 10 
Croston et al. 2008). 


5. Discussion 

Our XMM-Newton survey of galaxy clusters (Snowden et al. 2008), provides a large 
sample and high signal-to-noise data to study the scaling relationships between cluster mass 
and cluster X-ray temperature, gas mass, gas mass fraction, and the mass proxy Yx- Our 
fits are in good agreement with other observational and theoretical work, with one major 
caveat. We find a significantly larger scatter around the best-fit relationships than has been 
previously seen. 

In our sample, the scatter around the best-fit scaling relationships is 2-3 times the scat- 
ter found in most other observational and theoretical work. A recent study of the scaling 
relationship between the Compton y-parameter from SZ studies and cluster masses obtained 
from gravitational lensing also showed large scatter in the Mql — Tx and Mql ~ Y relation- 
ships (41% and 32%, respectively; Marrone et al. 2009). While other work has suggested 
that disturbed clusters will show more scatter than relaxed clusters (e.g. Kravtsov et al. 
2006), those authors used a visual classification of their simulated clusters, rather than a 
quantitative determination, making comparison difficult. Within our sample, we find no 
correlation between measures of the cluster dynamical state (ellipticity and asymmetry) and 
degree of scatter. 

Only one cluster property showed a strong correlation with degree of scatter, the combi- 
nation 0 e ff + pt which hydrostatic mass estimates, like those used here, depend on explicitly 
(e.g. Vikhlinin et al. 2006). Our data is the first to show this dependence due to the large 
sample size and range of cluster properties included. We note that there is no reason to 
suspect that true cluster masses would show such a dependence given the bias (and scatter) 
in hydrostatic mass determinations found in simulations where the true cluster masses are 
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known (see e.g, Nagai et al. 2007b). 

One question we must ask is how reliable are our determinations of the temperature 
profiles, whose wide range dominate the measurement of /? e // + (3 t . To check the reliability of 
the background modeling and spectral fitting procedure. Snowden et al. (2008) compared the 
profile of A 1795 with published Chandra , XMM-Nevuton, and Suzaku temperature profiles 
and found no significant difference. They also found reproducibility for three clusters with 
multiple observations. An initial comparison of our average cluster temperature profile is 
in good agreement with previous results both in overall shape and expected scatter (see 
Juett et al. 2009, in prep). The results of Leccardi & Molendi (2008) also suggest that 
cluster temperature profiles, while generally showing a falling profile at large radii, do show 
a range of profiles. Followup observations with Suzaku of the most unusual systems should 
be performed to confirm our results. 

Assuming our temperature profile range is indicative of the cluster population, we then 
need to ask how does this result affect cosmological studies using clusters. First, a larger 
scatter should be taken into consideration when discussing systematic errors in cosmological 
studies (see e.g. Vikhlinin et al. 2009). However, given other error sources, it is not clear 
that the mass scaling relationship scatter would be the dominant error contributor. 

It may be possible to correct for the scatter from theoretical studies of cluster properties 
but it is unclear if present models are consistent with our results. Nagai et al. (2007a) find 
little variation in the temperature profile at r 50 o in their simulations, although these are 
limited to their relaxed subsample. If a more thorough study of the simulations is not able 
to reproduce the observed cluster variation, that may point to some missing physics needed 
to better describe the conditions within clusters. 

Another result we would like to highlight is the comparison of our average gas mass 
fraction with previous results. The differences between our work and others (e.g., Vikhlinin 
et al. 2006; Allen et al. 2008) range from 5-20%. This is comparable with expected systematic 
differences between the instruments and analysis methods, but is significantly larger than the 
statistical errors typically quoted. In our opinion, a study of the expected systematics due to 
(1) instrumental calibration differences, and (2) data analysis methods must be performed. 
These issues are beyond the scope of this work, but we will address them in our future paper 
(Juett et al. 2009, in prep). 
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a contract with NASA. 
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Fig. 1. — Upper Left: Plot of X-ray spectral temperature, T x , and total cluster mass, M. 
Overplotted are the best-fit power-law relations using the BCES orthogonal slope estimator 
(red), the BCES Y|X slope estimator (blue), and a fixed slope of a = 1.5 (green). Upper 
Right: Plot of Y x and total cluster mass, M, with best-fit power-law relations overplotted. 
Color coding is the same as for the M — T x plot with a = 0.6. Lower Left: Plot of cluster 
gas mass, M g , and total cluster mass, M, and best-fit power-law relations overplotted. Color 
coding is the same as for the M — T x plot with a = 1.0. Lower Right: Plot of gas mass 
fraction, f g , and total cluster mass, M. Overplotted is the best-fit relation. 
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Fig. 2. — Plot of deviation (SM(Yx)/M; top panels) and scatter ({5M(Yx)/M} 2 ; bottom pan 
els ) around the best-fit M — Yx scaling relationship versus measures of the cluster dynamica 
state, ellipticity ( left panels) and asymmetry ( right panels). No correlation is found. 
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Fig. 3.— Plot of deviation around the best-fit mass scaling relations, SM(X)/M, vers 
(3 e ff + 0 t . The results for the scaling relationships for M — T x ( top panel ), M — M g ( midi 
panel) and M — Y x ( bottom panel) are given. All show a correlation between scatter ai 

fieff + / % ■ 




Table 1. Best-Fit Parameters of Powerlaw Fits 


Relation 8, 

Fit Method 5 

lOglO C/fgfl 

a 

Scatter 

M-T 

Y|X 

14.64±0.03 

1.67±0.16 

0.420 


Orth 

14.63±0.03 

1.89i0.20 

0.430 


Fix 

14.63±0.03 

1.5 

0.412 

M-Yx 

Y|X 

14.657±0.018 

0.61Y0.04 

0.218 


Orth 

14.657±0.018 

0.61±0.04 

0.218 


Fix 

14.646i0.018 

0.60 

0.212 

M-M, 

Y|X 

14.503±0.018 

0.93±0.05 

0.144 


Orth 

14.504i0.017 

0.93i0.05 

0.145 


Fix 

14.484i0.013 

1.0 

0.149 

/,- M 

Y|X 

0.134i0.005 

0.011i0.013 

0.104 


Orth 

0.134Y0.005 

0.011±0.013 

0.104 


a The relationships were of the form E(z) n M = C(X/X o)“ for the M — T, 
M — Yx, and M — M g fits with n = 1, 2/5, and 0, and X 0 = 5keV, 3 x 
10 14 M©keV and 4 x 10 13 M© for the T, Yx, and M g fits, respectively. For 
the f g — M fit, the relationship took the form f g = f g $ + o;log 10 (Af / 10 l0 M ©) . 

b We use the BCES fitting package and show the results from the Y|X and 
orthogonal fitting methods. In addition, we fit the data with a fixed slope, 
a , given by the expected self-similar relationships (a = 1.5, 1.0, and 0.6 for 
the T, Yx, and M g fits, respectively). 



